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ABSTRACT 

We model the broad-band emission from SNR RX J1713. 7-3946 including, for the first time, a 
consistent calculation of thermal X-ray emission together with non-thermal emission in a nonlinear 
diffusive shock acceleration (DSA) model. Our model tracks the evolution of the SNR including the 
plasma ionization state between the forward shock and the contact discontinuity. We use a plasma 
emissivity code to predict the thermal X-ray emission spectrum assuming the initially cold electrons 
are heated either by Coulomb collisions with the shock heated protons (the slowest possible heating), 
or come into instant equilibration with the protons. For either electron heating model, electrons reach 
> 10'' K rapidly and the X-ray line emission near 1 keV is more than 10 times as luminous as the 
underlying thermal bremsstrahlung continuum. Since recent Suzaku observations show no detectable 
line emission, this places strong constraints on the unshocked ambient medium density and on the 
relativistic electron to proton ratio. For the uniform circumstcllar medium (CSM) models we consider, 
the low densities and high relativistic electron to proton ratios required to match the Suzaku X-ray 
observations definitively rule out pion-decay as the emission process producing GeV-TeV photons. We 
show that leptonic models, where inverse- Compton scattering against the cosmic background radiation 
dominates the GeV-TeV emission, produce better fits to the broad-band thermal and non-thermal 
observations in a uniform CSM. 

Subject headings: acceleration of particles, shock waves, ISM: cosmic rays, ISM: supernova remnants, 
magnetic fields, turbulence 



1. INTRODUCTION 

The supernova remnant RX J1713. 7-3946 (G347.3-0.5) 
has been detected at photon energies ranging from radio 
to TeV 7-rays. The GeV-TeV detections in particular 
make this SNR an important test-bed for models of 
particle acceleration in astrophysical shocks, and a large 
number of fits to the data have been presented with an 
array of environmental and particle acceleration param- 
eters. Invariably, parameters are found that allow good 
fits to the non-thermal observations (or some sub-set 
of the observations). A critical question for cosmic-ray 
(CR) origin concerns the production of the GeV-TeV 7- 
rays. Are these 7-rays primarily from inverse-Compton 
(IC) emission from relativistic electrons, or pion-decay 
emission from the interaction of relativistic hadrons with 
the ambient medium? Models with good fits to the TeV 
emission with either inv erse-Compton or pion-decay 
have b een presented (e.g. , iPorter. Moskalenko fc Strong! 
2006t iBerezhko fc Volkl [2Q08I: iTanaka et al.l l2008r 
Morlino. A mato fc Blasi"2009': 'Zirakashv ili fc AharonianI 



and the thermal emission scales as rip. 

Until now, fits to the broad-band emission that incor- 
porate nonlinear diffusive shock acceleration (DSA) have 
not accurately accounted for the thermal X-ray emis- 
sion that might be present. We do this here for a SNR 
evolving in a uniform circumstellar medium (CSM) with 
no density enhancements as might occur with a pre-SN 
dense shell, nearby molecular cloud, etc. We find that the 
lack of observed thermal line emission eliminates pion- 
decay as the source of TeV emission in models with uni- 
form circumstellar media. 

The essential eleme nts of our CR-hydro-NEI mode l 
have been presented i n lEllison &: Cassam-Chena'J 



20091: lYam azaki. Kohri &: Katagiri 

200C~ 

and strong 



but conflicting claims for or against one or the other 
scenario, based on broa d-band continuum observations, 
have been made (e.g., iKatz &: Waxmaiil 120081 : iPlagal 
[2OOI IBerezhko fc"V6igi2009t) . We find that it is hard 
to discriminate on the basis of continuum emission 
alone, but that thermal X-ray line emission can easily 
differentiate between IC and pion-decay models because 
pion-decay requires a high proton number density, rip, 
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lEUison et aTl (|2007[) : iPatnaude. Ellison fc Slan 
and references therein. We couple a one-dimensional hy- 
drodynamic simulation of an evolving SNR with nonlin- 
ear diffusive shock acceleration. The ionization struc- 
ture, free electron number density, and electron tem- 
perature in the evolving interaction region between the 
forward shock (FS) and contact discontinuity (CD) 
are determined with a self-consistent treatment of the 
nonequilibrium ionization (NEI). We couple our com- 
puted nonequilibrium ionization fra ctions of heavy ele- 
ments to an updated version of the iRavmond fc SmithI 
(|1977[ ) plasma emissivity code to compute the thermal 
X-ray emission. 

Simultaneously, the shock accelerated, non-thermal 
electron and proton spectra are calculated, evolved, and 
used to determine the synchrotron, IC, non-thermal 
bremsstrahlung, and pion-decay emission from the SNR. 
We therefore obtain, for the first time, consistent thermal 
and non-thermal emission in an evolving SNR. 

2. MODEL 
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Any reasonably complete broad-band model of a SNR 
has a host of parameters. SNR RX J1713. 7-3946 is no 
exception and in this paper we do not present a full pa- 
rameter search. Instead, we concentrate on three essen- 
tial coupled components: (1) the SNR hydrodynamics, 
(2) nonlinear DSA, and (3) non-equilibrium ionization. 

Following the majority of work on SNR RX J1713.7- 
3946, we assume an age <snr — 1600 yr, and a distance 
Z?SNR — 1 kpc. Using the observed angular size, Dsnr 
implies a forward shock radius i?FS — 8.7 pc. While SNR 
RX J1713. 7-3946 is believed to be a core-collapse SN, we 
again follow the majority of work on this remnant and 
assume a uniform CSM with constant proton number 
density, rip, and constant unshocked magnetic field, Bq. 
We will present models where a prc-SN wind is assumed 
in future work. Besides Up and Bq, the following en- 
vironmental parameters are required to model the SNR 
evolution: the SN explosion energy, -Esn, the ejecta mass, 
Mej (we assume an exponential mass distribution for the 
ejecta), and the temperature of the unshocked CSM, Tq. 

We show models with two sets of parameters. In the 
"hadronic" model, the parameters are such that the GeV- 
TeV emission is dominated by pion-decay, while in the 
"leptonic" model, IC produces the GeV-TeV emission. In 
both cases, the parameters are chosen to simulta neously 
match the HESS TeV observations (fAharonian et al 



2007) and the Suzaku X-ray continuum (jTanaka et al 



2008|) . The hadronic and leptonic names refer to the 



particles, protons or electrons, mainly responsible for 
the GeV-TeV emission. As we show below, both models 
place the majority of the accelerated particle energy in 
protons, not electrons. 

We include an amplification factor, Samp, for the 
shocked magnetic field. In our simple ad hoc model 
of magnetic field amplification (MFA), the compressed 
magnetic field immediately behind the shock is in- 
creased by a factor i?amp- The amplified down- 
stream field is then ev olved in the downstream re- 
gion as described in lEllison et aTl (|2007[ ). For 
more self-consistent mode ls of MFA s ee, for exam- 
ple, IVL adimirov, Elfison fc Bvkovl (I2006D: ICaprioh et al.l 
(|2008l )r iVladimirov. Bvkov. fc EllisonI (j2009t) . 



To model the nonthermal radiation, we need additional 
parameters for nonlinear DSA.^ These are the accelera- 
tion efficiency, ^dsa (i-e., the instantaneous fraction of 
shock ram kinetic energy fiux placed in superthermal pro- 
tons), the relativistic electron to relativistic proton ratio, 
Ke-pf^ the maximum energy the protons obtain E™^^, and 
a factor, Ofcut, characterizing the shape of the turnover 
region around E™^^. We determine E™'^^ by limiting the 

^ The model of nonlinear DSA we use here is based on the semi- 
analytic model develo ped by IBlasi. Gabici fc Vannonil | |2005|) and 
lAmato fc Blasil l|2006t l. In our implementation, we fix the accel- 
eration efiiciency rather than the injection fraction as is done by 
Blasi and co-workers. While this difference may have important 
consequences during the early stages of the S NR evolution when 
the F S Mach number is extremely large (see IBerezhko fc EllisonI 
I1999I V it makes no significant difference to the integrated spectra 
at the later times we show here. 

* Note that K^p sots the post-shock relativistic electron density 
given the post-shock relativistic proton density. The post-shock 
thermal electron density, which determines the bremsstrahlung 
continuum and the X-ray line emission, is set by the densities and 
ionization states of the post-shock hydrogen and heavier elements. 
The model parameters K^p and rip are independent. 



acceleration when the acceleration time matches the SNR 
age or when the upstream diffusion length matches some 
fraction, /sk, of the shock radius, whichever comes first. ^ 
The factor acut smoothes the particle spectrum around 
£;max niimickiug the effects of particle escape (see, for ex- 
ample, IZirakashvili fc Ptuskinll2008l) . The above param- 
eters are fu lly defined in lEllison. Decourchelle fc Ballet! 
(f200l and lEllison fc Cassam-Chenaj (|2005l ). The ef- 
ficiency of DSA has been directly measured at the 
quasi-parallel Earth bow shoc k with ^dsa ^ 0.25 
(jEllison. Moebius fc PaschmannI Il990t ). Indirect ev- 
idence, based on particular models, suggests that 
the efficiency in some young SNRs, at least in 
some regions of the FS, can be 50% or more (e.g. 



Volk. Berezhko fc Ksenofontovl2003l:IWarren et ai]|2005l: 



Helder et al.ll2009D . 

For the thermal X-ray emission, we assume cosmic 
abundances and compare two extremes for heating the 
initially cold electrons. The slowest possible heating 
is from Coulomb collisions and the fastest is instant 
equilibration between electrons and protons, presum- 
ably produced by wave-particle interactions. For shock 
speeds above ~ 1000 km s~^, it has been suggested that 
electrons are heated very rapidly to kT ~ 0.3 keV by 
lower hybrid waves, after which continued heating to 
kT ~ 1 keV proceeds through Coulom b collisions (e.g., 
iGhavamian. Laming fc Rakowskil |2007() . Since, as we 
show below. Coulomb collisions alone rapidly heat the 
gas to ~ 0.3 keV, any difference between lower hybrid 
wave heating and Coulomb heating would only be impor- 
tant for UV and optical lines, so pure Coulomb models 
are appropriate for the X-ray emission. 

It is important to note that, in our CR-hydro-NEI 
model for the interaction region between the CD and FS, 
including X-ray line emission only requires two additional 
assumptions. One is the CSM elemental abundance and 
the other is the electron heating model. For Type la SNe, 
and a wide range of low-to-moderate mass core-collapse 
SNe, it is rea sonable to assunie solar abundances for 
the C SM (•e.g.. lChiosi fc Maedeil[l98l iKudritzki fc Pulsl 
|2000() . The two heating extremes we consider cover all 
likely possibilities. 



3. RESULTS 
For our leptonic model we assume n 



p 



0.05 cm 



Bo = 3mG, W = 0.25, /Vcp = 2x10-^ B., 



■'amp 



1, 

/sk — 0.1, and ttcut = 1- For the hadronic model. 
Hp = 0.2cm-3^ Bo = 2^G, ^dsa = 0.5, K^p = 7x10-"^, 
Samp = 5, /sk = 0.05, and acut = 1- In both mod- 
els, the values for i?sN and Moj arc varied with Up to 
obtain Rps ~ 8 — lOpc at ^snr = 1600 yr. Thus, for 
the leptonic model, Esn = 1 x 10^^ erg and A/oj = 3 Mq, 
while the hadronic model uses i?sN = 2 x 10''^ erg and 
Mej = 1.4 Mq.^ For a particular Up, other combinations 
of EsN and Moj giving Rps — 8 — 10 pc at 1600 yr yield 
similar results. In all cases, we assume Tq = lO^'K.^ 

^ The diffusion length in the FS precursor is determined assum- 
ing "Bohm diffusion," where a particle's mean free path is on the 
order of its gyroradius. 

^ The value M^j = 1.4 Mq is not meant to imply that we believe 
SNR RX J1713. 7-3946 originated from a Type la supernova. 

^ As long as To < 10*^ K, the unshocked temperature only weakly 
influences our results. 
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Fig. 1. — The top four panels show the free electron density, 
rie, the temperature, the ionization age, net, and the magnetic 
field in a parcel of gas first shocked at 200 yr. In panels (A) and 
(C), the solid curves are for Coulomb equilibration and the dashed 
curves (barely visible) are for instant equilibration. In panel (B), 
the dashed curve is the proton temperature and the solid curve 
the electron temperature assuming Coulomb equilibration. The 
dotted red curve in (B) shows the equal electron and proton tem- 
peratures assuming instant equilibration. In the bottom panel, the 
solid curves show the total emitted flux, at tsNR = 1600 yr, per 
arbitrary unit mass, in the band 1 — 2 keV from parcels of gas 
shocked at previous times. The dashed curves in the bottom panel 
show the corresponding flux from the bremsstrahlung continuum. 
The parameters are for our hadronic model with Up = 0.2 cm^^. 



At the end of the simulation, we obtain for the 
leptonic (hadronic) modeh the forward shock radius 
i?FS — 9.3 (8.8) pc; the forward shock speed Vps — 
3000 (2300) km s~^; the magnetic field immediately be- 
hind the FS i?2 — 10 (36) /iG; the overall FS compression 
ratio i?tot — 4.6 (5.6), the subshock compression ratio 
i?sub — 3.98(3.86), the fraction of SN explosion energy 
placed in CR ions ~ 0.13 (0.4), and the mass swept up 
by the FS ~6(19)Mo. 

In Fig. [T]we illustrate the properties of our CR-hydro- 
NEI model by following particular parcels of plasma. In 
the top four panels we show, for our hadronic model, the 
free electron number density, n,,, the electron and proton 
temperatures, the ionization parameter or age, net (t is 
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Fig. 2. — The top two panels show fits to the Suzaku RX J1713.7- 
3946 observations with our hadronic model for both Coulomb and 
instant temperature equilibration but ignoring the X-line emission. 
The blue (heavy wt.) solid curve is the synchrotron continuum, 
the black solid curve is the thermal emission (only lines above 
10~* MeV are included), and the dotted curve is the underlying 
bremsstrahlung continuum. The observed emission would be the 
sum (not shown in the top two panels) of the solid black and blue 
curves. The bottom panel shows the leptonic model (with Coulomb 
equilibration) where parameters have been chosen to be consistent 
with the Suzaku observations. For the hadronic model, the ra- 
diation intensity is multiplied by 0.95 to match the observations. 
For the leptonic model, a normalization factor of 0.2 is required to 
match the observations. We note that the Suzaku data have been 
adjusted for interstellar extinction so no extinction is applied to 
the model in this plot. 

the time since the parcel was shocked) , and the magnetic 
field in a parcel of plasma that is overtaken by the FS at 
200 yr. The red dotted curve in panel (B) gives the tem- 
perature assuming instant equilibration. Even though 
Te/Tp < 0.1 throughout the simulation for Coulomb 
equilibration, T^, approaches 10^ K (^ 850 eV) rapidly 
before leveling out. 

In the bottom panel of Figure [1] we plot the thermal 
X-ray emission between 1 and 2 keV, for both instant 
and Coulomb equilibration, at the end of the simulation 
for parcels of plasma shocked at previous times. The 
dashed curves are the continuum emission between 1 and 
2 keV and the total emission (solid curves), including 
lines, stands well above this regardless of the electron 
equilibration. As the left end of the bottom panel shows, 
at isNR — 1600 yr, plasma that was shocked > 200 years 
earlier is sufficiently ionized to produce a substantial flux 
in lines regardless of the electron equilibration. 
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In the top two panels of Figure [2] we compare 
our hadronic mode l to Suzaku observations of J1713 
(|Tanaka et al.1 12008[) for Coulomb (top panel) and in- 
stant equilibration (middle panel) . The Suzaku observa- 
tions have been adjusted for interstellar extinction and 
all model parameters are the same as in Fig. [1] For our 
hadronic model, we have chosen parameters that result 
in pion-decay dominating the GeV-TeV emission, i.e., Up 
must be above some limit and K^^ must be below some 
limit for this to be the case. Figure [5] makes it clear, 
however, that the X-ray line emission is much stronger 
in the hadronic model than can be accommodated by ob- 
servations. This is true for Coulomb equilibration even 
though the bremsstrahlung continuum remains well be- 
low the Suzaku observations. The only way to lower this 
emission relative to the synchrotron continuum would be 
to increase i^ep or to decrease rip to values that would 
then no longer reproduce the observed gamma-ray emis- 
sion. This is true regardless of the electron equilibration. 
Wc note that lowering in uniform CSM models re- 
quires lowering i?sN to maintain i?FS ~ 8 — 10 pc. 

We are unable to find any set of parameters that gives 
pion-decay dominating the TeV emission without pro- 
ducing emission lines around 1 keV that are inconsistent 
with the Suzaku observations. 

In the bottom panel of Figure [2] we show our Icptonic 
model where we have chosen parameters to be consis- 
tent with the smooth Suzaku observations. In addition 
to the parameters discussed already, we have arbitrar- 
ily adjusted the overall normalization of both models to 
match the observations. The hadronic model has been 
multiplied by 0.95 and the leptonic model by 0.2. Nor- 
malization values < 1 might correspond, observationally, 
to a partially complete shell morphology for the SNR, 
or possibly some reduction in the DSA injection and/or 
accelerati on efficiency over some f raction of the SNR sur- 
face fe.g.. lBerezhko fc V61kl[2008h . 

In Figure 3, we show our best fit hadronic and lep- 
tonic models, folded through the Suzaku XIS instrument 
response. * For both models, we simulated 20 ks obser- 
vations of the entire SNR with no background subtrac- 
tion, assuming a Galactic column density njj = 7.9 x 
10^^ cm~^. When co mpared to the Suza ku observations 
(cf.. Figs. 10 or 11 in lTanaka et al.ll2008D . it is clear that 
Suzaku would have detected lines as strong as those pro- 
duced in our hadronic model had they been present. 

In Fig. 21 we show broad-band fits to radio, Suzaku, 
preliminary Fermi- L AT, and HESS observations of RX 
J1713. 7-3946. The hadronic and leptonic models both 
produce reasonable fits if the thermal X-ray line emission 
is ignored. When the thermal X-rays are considered, the 
hadronic model is excluded. Only the cosmic microwave 
background is used to determine the IC emission. 

It is important to note in considering Figs.[5]and|3]tliat 
equally good fits to the continuum observations can be 
obtained with different parameter combinations. This, 
and the fact that the various models that have been 
applied to RX J1713. 7-3946 differ in details, accounts 
for the relatively small differences in parameters we ob- 
tain compared to those obtained by other modelers (e.g . , 
iBerezhko fc Volkl 120081: iMorlino. Amato fc Blasil l2009h . 

^ Response matrices are available at 

|http: / /heasarc. nasa.gov/docs/suzaku/prop_tools/xis_mat.htmil 




Energy [keV] 

Fig. 3.— Simulated Suzaku XIS spectra of RX J1713.3-3946. 
In the top panel, the best fit hadronic model is shown, with rip = 
0.2 cm^'^, while in the bottom panel, the best fit leptonic model 
is shown, with rip = 0.05 cm""^. In both panels, the blue curve is 
the contribution from the thermal X-ray emission, while the red 
curve is the contribution from synchrotron emission. The spectra 
correspond to a simulated 20 ks observation and are normalized to 
match the unabsorbed 1.0 - 10.0 keV fiux of 7.65 xlO"^" erg cm~^ 
s~^ found by Tanaka et al. (2009). In these simulated observations, 
we assume a Galactic njj = 7.9 X lO'^^ cm~^. 

However, consistency with the thermal X-ray line emis- 
sion forces the CSM density down and ifep up so no 
set of parameters can be found that result in pion-decay 
dominating the GeV-TeV emission. 

Characteristically of efficient DSA, the CR-hydro-NEI 
model produces an overall shock compression, i?tot > 4, 
and a subshock compression, i?sub < 4. Neverthe- 
less, even with 50% efficiency (5dsa = 0.5), i?sub re- 
mains large enough for electrons temperatures to be high 
enough for strong line production. 

The only factor we see that could lower the thermal 
emission substantially in a uniform CSM model, is the 
abundance. If the CSM is nearly devoid of heavy ele- 
ments, thermal line emission will be suppressed. Deple- 
tion onto dust will cut down C, Mg, Si and Fe, but it 
will not affect the O lines, which are the brightest in the 
model, or N or Ne. Furthermore, a substantial fraction 
of the dust is destroyed once net becomes a few times 
lO^^s-cm""^, so so me of the refractory elements would 
be liberated (e.g.. IWilliams et"al]|2006[ ). One does not 
expect really severe depletion in the low density uniform 
medium, but there could be significant dust in a red giant 
wind. 

It is also possible that the progenitor was a Wolf-Rayet 



5 



Hadronic 




-12 -9 -6-3 3 
Log^Q Energy [MeV] 



Fig. 4. — Broad- band fits to radio I IAcero et al.l 12009). Suzak u 
IITanaka et al.|[2003) . pre liminary Fermi- L AT lIFunk etaU [2009f) . 
and HESS observations IIAharonian et al.l 1200711 of RX J1713.7- 
3946. The top panel is our hadronic model and the bottom panel 
is our leptonic model. In both cases, the blue curve is synchrotron, 
the black is pion-decay, the red is IC, and the dotted is non-thermal 
bremsstrahlung. The dashed black curve is the sum including the 
X-ray line emission. As in Fig. [2] a normalization factor of 0.95 
(0.2) has been applied to the hadronic (leptonic) model. 



c 

V) 

o 



94 - 



92 



90 - 



88 



o 



1 1 1 1 1 1 


Solid: Hadronic _ 


/' 1 
e's // I 


Dashed: Leptonic. 






f 




-1 






-6 -3 


3 6 



Log.|Q Momentunn [m c] 



Fig. 5. — Phase-space distribution functions multiplied by for 
our hadronic (black solid curves) and leptonic (red dashed curves) 
models. These are integrated spectra at the end of the simulation 
in cgs units and represent the total material swept-up by the FS. 



(WR) star, and this could give anomalous abundances 
(e.g., lCrowtlieill2007[ ). Conversion of H to He reduces 
the number of electrons, weakening the line emission per 
unit mass by as much as a factor of two. However, WC 
and WO stars show much larger overabundances of O 
and Ne, which produce the strongest lines in the spectra, 
so the lines would be strongly enhanced. In WN stars, 
carbon and oxygen have been converted to nitrogen. The 
O lines would be weakened, and the N VH line at 500 
eV would be luminous but badly attenuated. The Ne 
IX and X lines at 922 and 1022 eV would then be the 
strongest in the spectrum at 0.5 to 1 times the strengths 
predicted. Thus, even in the case of a progenitor wind 
with anomalous abundances, we would still expect to see 
strong line emission in the swcpt-up CSM, and this would 
be present in the Suzaku observations. 

4. DISCUSSION AND CONCLUSIONS 

While several authors have proposed that emis- 
sion lines could be undetectable in J1713 because of 
low s hock temperature s or time-dependent ionization 
(e.g.. iDrurv et al.ll2009t iMorlino. Amato fc Blasil 120091: 
iBerezhko fc Volkl 12009( 1 we find that a SNR with prop- 
erties typically ascribed to J 1713, expanding in a uni- 
form CSM with solar elemental abundances, will produce 
strong X-ray emission lines when electron equilibration 
and non-equilibrium ionization are taken into account. 
This places constraints on the CSM density, n^, and on 
the relativistic electron to proton ratio, i^ep, to be con- 
sistent with Suzaku observations which show a smooth 
X-ray synchrotron continuum with no lines. 

While particular values of Up and Kep will depend 
somewhat on details of various DSA and SNR models, in 
any uniform CSM model the CSM must have a relatively 
low density and the electron to proton ratio of shock ac- 
celerated particles must be relatively high in order to pro- 
duce a satisfactory fit to the Suzaku data. Models where 
pion-decay produces the observed TeV emission require 
densities that are too high and values of Kep that are too 
low to be consistent with the Suzaku observations. We 
note that we have actually only computed a lower limit 
to the line emission since we have not included line emis- 
sion from the ejecta material heated by the reverse shock. 
If emission from a RS had been included, our conclusion 
that pion-decay is excluded could only be strengthened. 

Apart from minor differences, our fit to the broad- 
band spectrum (b ottom panel Fig. [H) is con s istent 
with others (e.g.. iPorter. Moskalenko fc Strong! 120061 : 
IMorlino. Amato fc Blasil l2009f ) where IC dominated the 
TeV emission. Our results differ substantially from the 
conjecture made by IDrurv et al.l (j2009D that the post- 
shock temperature can be reduced below X-ray emit- 
ting temperatures in strong shocks. The conclusions of 
IDrurv et all (|2009t ) are based on scaling arguments in 
the limit of extremely high sonic and Alfven Mach num- 
bers where the acceleration efficiency approaches 100%. 
In this case, the subshock m ay become weak en ough 
to limit heating to the values IDrurv et al.l (|2009D sug- 
ges t. However, Mach numbers as high as assumed in 
the IDrurv et al.l (|2009f ) scalings are not obtained for rea- 
sonable ambient magnetic fields and other parameters 
normally assigned to RX J1713. 7-3946. When nonlinear 
effects ar e fully taken into account for J 1 713 parameters 
(see also. IMorlino. Amato fc Blasil l200"9l ). and for accel- 
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eration efficiencies even modestly below 100%, the post- 
shock plasma (i.e., the proton com ponent) is heated more 
strongly than lDrurv et al.l (|2009[ ) suggest. 

Wc further emphasize that there is little freedom to 
reduce the thermal emission since we have calculated the 
NEI for the two electron heating extremes: Coulomb col- 
lisions and instant equilibration. For both extremes, and 
all cases in between, the shock heated plasma produces 
strong line emission. As Fig. [T] shows, it is not necessary 
for electrons to equilibrate with protons to become hot 
enough for line emission. For Coulomb collisions in our 
hadronic model, T^/Tp remains less than about 0.1 for 
> 1600 yr. 

Once it becomes clear that X-ray emission lines will 
be produced efficiently with a luminosity approximately 
oc rip, the intensities, /, of all the emission processes can 
be roughly scaled with the important parameters, Up, 
Kepi and the average downstream field B2 as: 

he oc nf oc nfK.p ; (1) 

Isyn « nfBf oc nfK^pBf ; (2) 

-^lino OC /brom OC 71^ ; (3) 

and 

Ipp oc Up^^Up oc -Faceup • (4) 

Here, the superscript "rel" indicates the number density 
of relativistic particles capable of producing the observed 
radiation. The factor i^acc is some function of the DSA 
efficiency, i.e., the fraction of ambient protons turned into 
relativistic protons capable of producing GcV-TeV emis- 
sion (n™' oc Faccnp). We also assume that the relativistic 
protons producing pion-decay are drawn from the same 
population as the target protons. The expression for /gyn 
assumes the underlying electron spectrum is a power law, 
dN/dE oc E^", with cr = 2a - 1. 

If the TeV emission is from pion-decay, then the ratio 
Jsyn//pp is fixed by the observations and 

^syn/^pp oc KcpB^/Up = G , (5) 

where G is some constant determined by either the radio 
or X-ray synchrotron observations. If G is set by radio 
observations, radiation losses don't play a role. To hide 
the X-ray lines, we need to increase 

/syn//ii„c cc nfK,pBf/nl = nfG/np cx FaccG . (6) 

Thus, the only parameter that can change the relative in- 
tensity /gyn/^iino is thc DSA efficiency. The X-ray line to 
synchrotron continuum ratio can be changed by changing 
the magnetic field, but the absolute ratio of X-ray lines 
to gamma rays is basically fixed in the hadronic scenario. 
From equations 3 and 4, 

hnc/Ipp OC 1/Facc , (7) 

and for the hadronic model (top panel in Fig.|4]), /iinc/-^pp 
is more than an order of magnitude too large compared 
to observations to be accommodated. Changing Face 
and/or B cannot hide the lines if Up is too large. 

Of course, the situation is more complicated for several 
reasons. (1) The line emission depends importantly on 
the SNR evolution (i.e., the ionization age; Fig. [Ij and 
the CSM composition. (2) Thc factor Face depends on 



the shock dynamics, the magnetic field, and uncertain 
details of NL DSA. Furthermore, since radiation losses 
are important for relativistic electrons but not relativis- 
tic protons, n^'/rip*^' 7^ K^p at high energies and K^p 
depends on B for X-ray synchrotron emission. This will 
change the /syn/^iinc scaling at X-ray energies but will 
not change the relative intensities of radio vs. pion-decay 
emission. (3) The detailed fits to the Suzaku, Fermi-LAT, 
and HESS data depend critically on the shape of the 
underlying electron and proton spectra in thc turnover 
region, and on the SNR magnetic field morphology. De- 
spite these complications, Eqs. ^ and (O, must largely 
control thc overall scaling. 

For our hadronic model shown in Figs. and |31 we have 
chosen particular values of rip, K^p, and i?amp to match 
the shape and relative normalization of the radio. X-ray, 
and TeV observations. For the acceleration efficiency, we 
have used ^dsa = 0.5, i.e., 50% of the instantaneous FS 
ram kinetic energy flux is put into relativistic protons. 
While there is little indication that larger values of ^dsa 
occur in SNRs, wc have explored ^'dsa > 0.7 and find 
a poorer match to the broadband observations and no 
improvement in the hadronic fit to the X-ray lines. One 
reason for this is that, in nonlinear DSA, an increase in 
acceleration efficiency must be accompanied by an in- 
crease in the overall shock compression ratio, i?tot- This 
translates to an increase in the downstream plasma den- 
sity, a decrease in the electron temperature equilibration 
time, and stronger X-ray line production. Furthermore, 
increasing the acceleration efficiency also increases B2 
due to compression and possibly more by MFA. Because 
of changes in i?2, increases in ^dsa are constrained by 
equation [S] 

On thc other hand, it is easy to show that lowering 
^DSA below 0.1 is also inconsistent with the broad-band 
observations. 

We emphasize again that the modeled shape of the 
high-energy turnover is both critical and uncertain. For 
IC and synchrotron, the shape depends on the compe- 
tition between acceleration and radiation loss timescales 
in the acceleration region. The turnover will be further 
modified by radiation losses as the electrons evolve be- 
hind the shock and by diffusion of the high-energy elec- 
trons into regions of different density and magnetic field. 
In fact, high-energy electrons might diffuse away from re- 
gions of high magnetic field, reducing their synchrotron 
emission while they still emit IC, For pion-decay, the 
turnover depends on the maximum energy the FS can 
produce which depends on the self-generated diffusion of 
thc highest energy, escaping particles. 

Since these effects are yet to be described precisely, 
all existing SNR models, including ours, make arbitrary 
approximations that importantly influence the turnover 
shape. The fit to the shape of the Fermi-LAT and 
HESS observations is determined largely by /sk, ctcut, 
and -Bamp- The detailed fit to the shape of the Suzaku 
observations depends largely on ctcut and -Bamp- 

Other effects may be important as well. The 
synchrotron spectrum might be hardened in 
thc turnov er region by stochastic e ffects , as de- 
scribed in iBvkov. Uvarov. fc Elliso^ (|2008D . Fur- 
thermore, as suggested by sev eral authors (e.g., 
iPorter. Moskalenko fc StronsJ |2006( ) , a photon source 
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in addition to the cosmic microwave background might 
improve the IC fit to the highest energy HESS points. 

In contrast to the shapes of the radiation spectra near 
their maximum energies, the relative normahzations of 
synchrotron, IC, pion-decay, and thermal X-ray emission 
are less uncertain because they depend more concretely 
on basic parameters. We believe that none of the ap- 
proximations in our CR-hydro-NEI model are significant 
enough to change our basic conclusion: the constraints 
on ambient density and K^p from thermal X-ray emission 
rule out pion-decay as the mechanism producing TeV 
emission in models with a uniform CSM. 

The fact that electrons are likely producing the high- 
est energy photons observed from RX J 1713. 7-3946, does 
not lead us to suggest that protons arc absent or less en- 
ergetic. Our Icptonic model assumes that at any instant 
25% of the shock ram kinetic energy flux goes into rela- 
tivistic protons while less than 1% goes into relativistic 
electrons. Electrons are observed simply because they 
radiate more efficiently than protons in low density me- 
dia. As Fig. [S] shows, the maximum proton energy is 
similar in our two models, i.e., E^'^^ ~ 10"^'' eV. The 
increase in E^'^^ from the higher shocked magnetic field 
in the hadronic model (e.g., B2 — 36 at the end of 
the simulation vs. B2 — 10 /iG for the leptonic model) 
is partially offset by the smaller /sk factor (/sk = 0.05 
for the hadronic model while /sk = 0.1 for the leptonic 
model). The electron maximum energy is about a fac- 
tor of 10 higher in the leptonic model, and the shapes 
of the electron spectra are different, due to the effects of 
radiation losses. 

One result of our leptonic model, which integrates 
emission over the entire remnant, that may conflict with 
observations is the low shocked magnetic field. A low 
B2 favors the leptonic model and we obtain B2 — 10 /iG 
for the parameters used here. Much higher estimates 
for B2 have been obtained for thin X-ray filaments 
where the sharp X-ray edges and/or rapid time varia- 
tions are attributed t o strong synchrotron losses (e.g., 
lUchivama et al.ll2007D . Our uniform CSM assumption 
cannot descibe filaments and it is possible that more com- 



plicated, multi-component models could account for this. 
For example, if the synchrotron emission originates from 
a smaller region than the IC emission (due, for exam- 
ple, to a strong but compact postshock magnetic field), 
then a larger field strength would be possible for a given 
invcrsc-Compton flux. 

We have been careful to emphasize that we only con- 
sider a uniform CSM in this paper. While pion-decay is 
eliminated in this simplest case, SNR RX J1713. 7-3946 is 
certainly more complex. As fermi-L^T" observations im- 
prove with time, the shape of the combined Fermi-LAT 
and HESS observations may indicate that the GcV-TeV 
emission is, in fact, hadronic in origin. This will require 
some multi-component model where relativistic protons 
interact with a high density target but care must still 
be taken to avoid inconsistency with the Suzaku obser- 
vations. If the FS runs into a high density shell, strong 
X-ray lines are likely to be produced along with the en- 
hanced pion-decay emission. If the highest energy pro- 
tons escape upstream from the FS and impact a high den- 
sity medium before the material is shock heated, pion- 
decay emission may be strong without strong accom- 
pany ing X-ray line emission (e.g.- lLee. Kamae fc EllisonI 
[2OO8I) . 
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